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Abstract: The international Muon Ionization Cooling Experiment (MICE), which is under con- 
struction at the Rutherford Appleton Laboratory (RAL), will demonstrate the principle of ionization 
cooling as a technique for the reduction of the phase-space volume occupied by a muon beam. Ion- 
ization cooling channels are required for the Neutrino Factory and the Muon Collider. MICE will 
evaluate in detail the performance of a single lattice cell of the Feasibility Study 2 cooling channel. 
The MICE Muon Beam has been constructed at the ISIS synchrotron at RAL, and in MICE Step I, 
it has been characterized using the MICE beam-instrumentation system. In this paper, the MICE 
Muon Beam and beam- line instrumentation are described. The muon rate is presented as a function 
of the beam loss generated by the MICE target dipping into the ISIS proton beam. For a 1 V signal 
from the ISIS beam-loss monitors downstream of our target we obtain a 30 KHz instantaneous 
muon rate, with a neglible pion contamination in the beam. 
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1. Introduction 

Muon storage rings have been proposed as sources of intense, high-energy neutrino beams flU, 0] 
and as high energy lepton-antilepton colliders [0, H]- In each of these facilities, the volume of phase 
space occupied by the muon beam must be reduced (cooled) before the beam is accelerated and 
stored. The muon lifetime is so short that unacceptable decay losses will occur if conventional 
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cooling techniques are employed [0]. Ionization cooling is fast, as the cooling is generated by ion- 
ization energy loss as the muon beam passes through an "absorber" material, the lost longitudinal 
momentum being restored in accelerating cavities. Ionization cooling has therefore been adopted 
in the proposed design for both the Neutrino Factory and the Muon Collider [0]. 

At the Neutrino Factory [0], intense beams of high energy neutrinos are produced from the 
decay of muons circulating in a storage ring. Long straight sections in the storage ring direct 
neutrino beams to one or more distant detectors. The energy spectrum and flavour composition 
of the beam are known precisely, as compared to those of conventional neutrino beams ED- A 
number of conceptual designs for the Neutrino Factory have been proposed [0, 0, HO, [H, Ol each 
of which exploits an ionization-cooling channel to increase significantly the muon flux delivered to 
the storage ring. The baseline design for the facility being developed by the International Design 
Study for the Neutrino Factory (IDS-NF) collaboration delivers 10^^ muon decays per year and has 
been shown to have a discovery reach that is significantly better than that of realistic alternatives. 
The cooling channel specified in the IDS-NF baseline delivers an increase of a factor 2.4 in the 
stored-muon flux and is essential for the design performance of the facility to be delivered. 

The Muon Collider flTlll offers an attractive route to multi-TeV lepton-antilepton collisions. 
Since the muon mass is 200 times that of the electron it is possible to deliver very high-energy 
beams using circular accelerators that can be designed to fit within the area of existing particle- 
physics laboratory sites. In addition, the large muon mass implies a reduction in the rate of beam- 
strahlung by a factor of ^ 10"^ over an e^e~ collider of the same centre-of-mass energy. As a 
consequence, the annihilation-energy distribution is much narrower at a Muon Collider than at an 
e^e~ collider of the same energy. 

The \i^\i~ annihilation cross section falls rapidly as the centre-of-mass energy in the colli- 
sion increases, therefore high luminosity (> 10^"^ cm~^s~^) is critical to the success of the Muon 
Collider. To achieve the required luminosity, the muon-beam phase space must be cooled in all 
six phase- space dimensions: a reduction in the six-dimensional emittance by a factor ^10^ being 
required. By contrast, the cooling requirements of the Neutrino Factory, which requires cooling of 
only the four-dimensional transverse phase space, are relatively modest. 

Ionization cooling is accomplished by passing the muon beam through a low-Z material (the 
"absorber") in which it loses energy through ionization. The energy loss results in a reduction 
in the longitudinal and the transverse components of momentum. The lost energy is restored by 
accelerating the beam such that the longitudinal component of momentum is increased while the 
transverse components remain unchanged. The net effect of these two operations is to reduce the 
divergence of the beam, thereby reducing the volume of transverse phase space it occupies. Beam 
transport through the absorber and the accelerating structures is achieved using a solenoid- focusing 
lattice. 

The rate of change of normalised emittance due to ionization cooling in a medium of thickness 
X may be described as [0]: 

dE^ ^ £n /dE\ ^ j8,(0.014GeV)^ ^ 



dX p^E^\dx/ ip^E^m^Xo 

where £n is the normalised transverse (four-dimensional) emittance of the beam, Xq is the radiation 
length of the medium, jS^ is the betatron function, E"^ and m^ the energy and mass of the muons and 
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j8 = pc/E their velocity. The first (negative) term on the right hand side of equation |0| describes 
a reduction of emittance per unit length, i.e. cooling. The second (positive) term describes the 
heating effect of multiple scattering. The two effects reach an equilibrium when the emittance of 
the beam is: 



ft (0.014 GeV)2 I dE 



g.^- oo.,. .. {—.) ' (1-2) 



The ideal cooling channel will produce the lowest equilibrium emittance. This is obtained when ft 
is minimised and Xq • ( ^ ) is maximised. To minimise ft requires strong focusing at the absorber 
while hydrogen offers the largest value of Xq • (^). 

While the principle of ionization cooling is readily described using equations |0| and [T^ , the 
construction of an ionization-cooling cell requires significant engineering. In a practical ionization- 
cooling cell, the energy loss in the absorber must be replaced by means of RF cavities: momentum 
must be restored so that beam transport through subsequent cells is unaffected. 

The international Muon Ionization Cooling Experiment (MICE) OH, [T51] is being constructed 
at the Rutherford Appleton Laboratory (RAL). The experiment will consist of the full cell of an 
ionization-cooling lattice and the instrumentation necessary to measure the emittance of a muon 
beam before it enters and after it leaves the cell. In this way, the MICE collaboration will measure 
the cooling performance of the lattice cell in a variety of modes of operation and over a range of 
momentum and emittance. The results of the experiment will allow the cooling channels of the 
Neutrino Factory and Muon Collider to be optimised. 

2. The MICE concept 

A schematic diagram of the MICE experiment is shown in figure |T|. The MICE cooling channel, 
which is based on a single lattice cell of the cooling channel described in [0], comprises three 
volumes of ^ 20 litre of liquid hydrogen and two linac modules each consisting of four 201 MHz 
cavities, with gradients of ^ 10 MV/m. Beam transport is achieved by means of a series of super- 
conducting solenoids: the "focus coils" focus the beam into the liquid-hydrogen absorbers, while 
a "coupling coil" surrounds each of the linac modules. 

Detailed Monte Carlo simulations of the passage of muons through the cooling cell have been 
carried out to estimate its performance. The results indicate that a reduction in the normalised 
emittance of 10% is expected for a muon beam entering the cell with a momentum of 200MeV/c 
and an emittance 8^ ^6.271 nmirad. The instrumentation upstream and downstream of the cooling 
cell is required to measure this change in emittance, Ae^^, with a relative precision Ae^/Sn =1%; 
i.e., measurements of £n upstream and downstream of the cooling cell with an absolute precision 
of 0.1% are required, to allow the extrapolation to a full cooling channel. Conventional emittance 
measurement techniques based on beam-profile monitors do not have the required precision. 

While the muon-beam intensity in the ionization-cooling channel at the Neutrino Factory and 
the Muon Collider is in excess of lO^'^/i^/s, the phase-space density is always low enough for 
space-charge forces to be neglected. This fact has been exploited to devise a single-particle exper- 
iment in which each muon is measured using state-of-the-art particle detectors and the bunched 
muon-beam is reconstructed offline. The instrumentation upstream of the MICE cooling cell 
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Figure 1. Top panel: view of the MICE experimental set-up; the cooling channel, with its three liquid 
hydrogen absorbers and two RF cavity modules, is sandwiched between two identical trackers. The sequence 
of solenoids defining the MICE optics is also visible. The muon beam is coming from the left. Bottom panel: 
3D cut-away engineering drawing of MICE, starting at first Spectrometer Solenoid. Beam instrumentation 
is not shown. 
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includes a particle-identification (PID) system (scintillator time-of-flight hodoscopes TOFO and 
TOFl and threshold Cherenkov counters Ckova and Ckovb) that allows a pure muon beam to be 
selected. Downstream of the cooling channel, a final hodoscope (TOF2) and a calorimeter system 
allow muon decays to be identified and rejected. The calorimeter system is composed of a lead- 
scintillator section (KL), similar to the KLOE flT^ design but with thinner lead foils, followed by a 
fully active scintillator detector (the electron-muon ranger, EMR) in which the muons are brought 
to rest. Charged-particle tracking in MICE is provided by two solenoidal spectrometers. The MICE 
instrumentation must be sufficiently robust to perform efficiently in the presence of background in- 
duced by X-rays produced in the RE cavities. Eor a full description of the experiment see flT511 . 

The position, momentum and energy of each muon is measured before and after the cooling 
cell. Eor each particle x, y, t, x' = dx/dz = Px/Pz^ / = dy/dz = Py/Pz and t' = dt/dz = E/p^ are 
measured; pt is the fl^ component of momentum, E is the energy and the MICE coordinate system is 
defined such that the z axis is parallel to the nominal beam axis, the 3; axis points vertically upward 
and the x axis completes a right-handed coordinate system. The time, t, corresponds to the moment 
the particle crosses a reference surface within the tracking volume. The input and output normalised 
emittances, 8^^ and 8^^, of a beam ofN muons can be determined from the measurements of the 
phase-space coordinates of single muons in the up- and down-stream trackers. Because they are 
measured for the same ensemble of muons, Ca^. and 8^^ will be correlated with the effect of reducing 
the statistical error on the fractional change in emittance, / = (Ca^. — 8nJ/8n, below \/\/N. The 
only stochastic process is multiple scattering in the absorbers and it can be shown that flTTII 

"'-IS «•'> 

where p^ is the mean square transverse momentum of the beam at the absorbers and p] is the mean 
square transverse momentum kick due to multiple scattering. It has been demonstrated with Monte 
Carlo studies flT^ that for beams with Ca^ > Stt mm • rad a sample of 3 x 10^ muons is sufficient to 
reduce the statistical error on / to less than one percent. 

The MICE Muon Beam is required to deliver beams with a pion contamination below 10% and 
a central momentum {p^) in the range 140 < p^ < 240MeV/c. The beam line must also deliver a 
transverse emittance in the range 3 < 8a^ < 10 tt nrni • rad. The full range of emittance is required 
over the full range of momentum. A tungsten or brass "diffuser" of variable thickness is placed 
at the entrance to the upstream spectrometer solenoid to generate the divergence necessary for the 
required range of emittance. 

At a rate of ^ 30/i/s a sample of 10^ muons may be accumulated in one hour, giving the 
possibility to record all the needed momentum/emittance beam settings in the available running 
time. 

The MICE program was conceived to be carried out in six "Steps", summarised in table [T| and 
figure 0. A thorough description of all MICE Steps is beyond the scope of this paper. However, each 
Step requires that the MICE Muon Beam meets the specifications outlined above. Step I, in which 
the muon beam was characterized using the time-of-flight (TOE), Cherenkov and KL detectors, 
was completed in August 2010. The EMR calorimeter flT9l1 . that will complete the downstream 
calorimetry, will follow soon. 
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Table 1. The six "Steps" in which it was conceived to carry out the MICE programme. The physics pro- 
gramme of steps II to III will be carried out within the STEP IV configuration. 



Step 



I 
II 

III 

IV 
V 

VI 



Description 



Beam characterisation and detector calibration. 
Upstream spectrometer and diffuser. Incoming beam emittance 
can be changed and measured precisely. 

Downstream spectrometer: study of systematics in emittance 
measurement in the two spectrometers. Solid absorber with mea- 
surement of energy loss and multiple scattering correlations. 
Introduction of the first liquid hydrogen absorber and focus-coil 
module. 

First linac section and second hydrogen absorber/focus-coil 
module. First measurement of cooling with momentum recov- 
ery. 

Second linac module and third liquid hydrogen absorber/focus- 
coil module. Measurement of emittance before and after the 
cooling cell, in various optics configurations, with momentum 
recovery. 







STEP I 




STEPV 






I 



STEP VI 



Figure 2. Schematic diagrams of the six "Steps" in which the MICE programme was conceived. 
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Figure 3. Top view of the MICE beam line with its instrumentation, as used in Step I. 



3. MICE Muon Beam 

3.1 Overview 

The design of the MICE Muon Beam is similar to that of the jUE4 beam line at the Paul Scherrer 
Institute (PSI) in Switzerland. The MICE Muon Beam may be divided, conceptually, into three 
sections (see figure ^. In the upstream section, pions produced by the ISIS proton beam striking 
a titanium target are captured using a quadrupole triplet (Ql-3). The pions are transported to a 
bending magnet (Dl) which directs pions of particular momentum into the decay solenoid (DS). 
The decay solenoid, the second section of the beam line, causes the pions to spiral around the 
nominal beam axis, increasing the path length. This effect together with the solenoid focusing 
increases the number of muons captured between Dl and the second dipole magnet (D2). D2, the 
first element of the downstream section, directs muons of a particular momentum into a quadrupole 
channel (Q4-6 and Q7-9) that transports the muon beam to the MICE experiment. 

The composition and momentum spectrum of the beam delivered to MICE is determined by 
the interplay between the two bending magnets Dl and D2. This is illustrated in figure 0. With 
Dl set to select a particular pion momentum, poi (horizontal axis), the kinematics of pion decay 
result in muons with momentum in the range Pj^i mm < P/i < PfimsLx- The kinematic limits are 
indicated in figure g. For pm ^ 200MeV/c, p^imax ^ Pdi- When D2 is set to select particles with 
momentum poi ^ Piimin, backward-going muons in the pion rest frame are selected. The muon- 
beam purity, P^, is defined as P^ = N^/Nj, where N^ is the number of muons in the beam and 
Nj is the total number of particles. A muon beam of sufficient purity can be selected by setting 
Pd2 = Piimin^ Pdi/2. This can be seen in the right-hand panel of figure ^ where the momentum 
spectra of particles emerging from the decay solenoid are plotted in a simulation in which Dl 
was set to select a central momentum of 400 MeV/c. The figure shows that setting D2 to select a 
momentum bite of 10% around pd2 = 236 MeV/c delivers a muon beam with a pion contamination 
(defined as Nji/Nj, where A^;^ is the number of pions in the beam) less than 2%. In practice, a muon 
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Figure 4. Working principle of the MICE beam line. Left-hand panel: the red and blue lines are the 
kinematic limits of the spectrum for muons produced in pion decays. By tuning D2 to the backward going 
muon peak an almost pion-free sample is produced. Right-hand panel: simulation showing pion and muon 
spectra at the end of the decay solenoid. Only high momentum pions survive. The green band shows the 
acceptance of D2, when tuned to the backward-going muon peak. 



beam with a purity in excess of 95% is obtained by setting pm — ^poi- This case is referred to as 
"71 -^ jU mode". A momentum- selected beam containing pions, muons, protons and electrons (or 
positrons) can be obtained by setting p^i — Poi- Such a setting is referred to as "calibration mode" 
and is used for the calibration of detectors. 



The muon beams delivered to MICE are characterized by their normalised emittance, Ca^, and 
momentum, p^. Different beam settings are identified using the notation {en^Pz), where £n is 
measured in k mm • rad and p^ in MeV/c. The "nominal" values of £n and p^ are defined such that 
the nominal value of 8^ is evaluated in the upstream spectrometer solenoid and the nominal value 
of Pz is evaluated at the centre of the central liquid-hydrogen absorber in the Step VI configuration. 
The baseline optics deliver a beam characterized by {En^Pz) = (6,200). The (En^Pz) space is 
discretised into a matrix of nine elements (see table ^. Each cell of this emittance-momentum 
matrix is associated with a specific downstream beam-line optics. Together the cells span the 
i^NiPz) space required to serve the entire MICE programme. 
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Figure 5. Left panel: MICE target drive schematics. Right panel: The MICE target installed in ISIS. 

3.2 Beam-line elements and instrumentation 

3.2.1 The MICE target 

Protons are injected into ISIS with a kinetic energy of 70MeV and are accelerated over a period 
of 10 ms up to 800 MeV at a repetition rate of 50 Hz ^ The MICE target has been designed to dip 
into the beam shortly before extraction at a rate of ^ 1 Hz, i.e. sampling one in 50 of the ISIS beam 
pulses. 

The MICE target is a hollow titanium cylinder with an outer diameter of 6 mm and an inner 
bore of 4.6 mm. The target dips vertically into the beam, intercepting it over the last 3 ms of the 
acceleration cycle (^ 4000 turns). Over this period, the target samples proton-beam energies in the 
range 617-800 MeV. The beam envelope shrinks as the beam is accelerated. 

To insert the target into the shrinking beam envelope requires an acceleration of ^ 80 g. The 
target-insertion mechanism consists of a linear electromagnetic motor in which permanent magnets 
mounted on the target shaft are accelerated by a series of 24 coils contained within the stator. An 
optical position- sensing system with laser quadrature readout is used to monitor the position of 
the shaft and to control the drive current to the 24 coils (see left-hand panel of figure ||). A paper 
describing the design and operation of the MICE target is in preparation fEDII . 

3.2.2 The MICE luminosity monitoring 

There are 39 beam loss monitors (BLM) fOTll placed around the inner radius of ISIS. These are 
argon-filled cylindrical ionization chambers, 3 m long, 16 mm in diameter and situated 2-3 m from 
the beam axis parallel to the beam. The BLMs detect particles produced by protons that are lost 
from the beam envelope and interact with accelerator components. The ionization chambers that 
form the BLMs produce a voltage signal proportional to the particle rate. This signal is integrated 
over some time interval (typically between 1 and 3 ms) to give an overall measure of beam loss 
(in units of V- ms). The ISIS synchrotron lattice has 10 "super periods", numbered from to 9, 



Upon reaching 800 Mev energy the protons are extracted to the ISIS neutron-production targets. 
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Luminosity monitor box (200 x 300 mm) 




85 mm 



PMT3 



75 mm 150 mm 



Figure 6. Left panel: Luminosity Monitor design schematic (particles from the beam travel from right to 
left). Right panel: photograph of the detector. 

with the MICE target situated in super period 7 (SP7). The signal from the four BLMs in SP7 is a 
measure of beam loss produced by the MICE target. 

The MICE Luminosity Monitor (LM) provides a measurement of the particle rate close to the 
MICE target independent of the ISIS Beam Loss Monitors. The LM is located inside the ISIS 
synchrotron ring, 10 m from the target (see figures ^ and ^. A line drawn from the target to the 
LM makes an angle of 25° with respect to the direction of the proton beam at the target. The LM 
therefore samples particles at an angle similar to that at which particles are captured by Ql-3. The 
rate recorded by LM may therefore be used to validate the simulation of the upstream beam line 
and to normalise the rates measured in the beam line instrumentation to the number of protons on 
target (POT). 

The LM consists of two pairs of small scintillators placed either side of a polyethylene block 
that absorbs protons with momentum below ^ 500 MeV/c and pions with momentum below ^ 
150MeV/c (see figure ^. The scintillators are read out by Hamamatsu H5783P low noise PMTs 
which have a rise time of 0.8 ns and a gain of 10^ fE2l1 . Upstream of the polyethylene block, PMTs 
1 and 2 are bonded to scintillator blocks of dimension 20 x 20 x 5 mm^ ; downstream, PMTs 3 and 
4 are attached to 30 x 30 x 5 mm^ scintillator blocks. Neutrons may be detected if scattering inside 
the polyethylene block releases a proton. The PMT signals are fed to LeCroy 623B NIM octal 
discriminators, the three outputs of which are used to form three coincidence signals using a LeCroy 
622 NIM quad coincidence unit fE311 . Coincidences are recorded within the experimental trigger 
gate during each spill: PMTs 1 and 2 (LM12), 3 and 4 (LM34) and all four PMTs (LM1234) ^. 

The quad LM1234 coincidence rate is plotted as a function of beam loss in figure 0. The data 
show an approximately linear relation between LM counts and beam loss in SP7 between 0.5 Vms 
and 3Vms. Due to beam losses not caused by the MICE target, the relation below 0.5 Vms 
becomes not linear. Above 3 Vms saturation effects become significant. The rate recorded by 



^In the following, we use the term "spill" for the burst of particles resulting from one target dip, even though the 
primary ISIS proton beam is not extracted in producing the MICE beam. 
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Figure 7. Luminosity monitor data showing a direct relationship between luminosity and beam loss on all 
three scalar output channels. Graphs show counts for coincidences of signals for PMTs 1,2,3 and 4. The 
straight lines are a linear fit to the data between 0.5 V- ms and 3 V- ms. 



the Luminosity Monitor has been used to convert integrated beam loss to protons on target (POT) 

3.2.3 The MICE Muon Beam magnets 

The quadrupole magnets that form the upstream pion-capture section of the upstream beam line 
are Type-IV quadrupoles recovered from the NIMROD accelerator which operated at RAL during 
the 1960s and 1970s fESII . Ql-3 have a circular aperture at the pole tips of 203 mm. The operating 
field gradients, shown in table §, are always significantly smaller than the maximum gradient of 
10 T/m. The power to each quadrupole is provided by a 200 A, 30 V power supply. An aluminium 
vacuum pipe passes through Ql-3 and extends into a vacuum box inside Dl. 

The quadrupoles that form the two downstream triplets (Q4-6 and Q7-9) are each Type QC 
quadrupoles from DESY fESH . The downstream quadrupoles have cross-shaped apertures, with a 
circular aperture of 352 mm diameter at the pole tips and a maximum gradient of 8 T/m. They are 
each excited by a 400 A, 70 V power supply. 

The two dipoles (Dl, D2) are large, 17 1, rectangular NIMROD Type-1 dipoles [E3]. The 
nominal apertures for these magnets are rectangular with tapered pole pieces that provide a 6 inch 
vertical aperture and a 20 inch horizontal aperture. Dl is operated at a field close to saturation of 
the iron yoke. The power supply for this magnet is limited for operation at a maximum current of 
440 A at 240 V. D2 is operated at a more modest field and is excited with a power supply delivering 
200Aatl00V. 

The large superconducting decay solenoid (DS) was supplied by the Paul Scherrer Institute 
(PSI) in Switzerland, where it had been in use at the /iE4 beam line from 1974 to 2004 [ETD. 
The coils which make up the DS have an open inner radius of 57.5 mm (outer radius 65 mm) and a 
length of 5 m. The conductor is a Cu/NbTi mixture with proportions Cu:NbTi of 3.5: 1 . The magnet 
has a stored energy of 1.5 MJ. The DS can operate at a maximum current of 1000 A, with a current 
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Table 2. Beam-line magnet physical properties. 



Beam line Magnet 


Physical length (mm) 


Aperture (mm) 


Typical field (limit) 


Quad Type-IV 


1143 


203 (/) 


1.6(10)T/m 


Quad Type QC 


1170 


352 (/) 


2 (8) T/m 


DipoleType-16" 


1400 


508 (H)x 152 (V) 


1.47 (1.5) T 


PSI Decay Solenoid 


5000 


115 (/) 


5 (5.7) T 



density of 220 A mm ^ . The nominal operating field in MICE is 5 T, corresponding to a current of 
870 A. 

A summary of the position and dimensions of all the components along the MICE beam line 
is presented in table 0. A photograph of the upstream section of the beam-line is shown in figure § 
and of the downstream beam line in figure 0. 

Table 3. The MICE beam-Hne elements and detectors for Step I. 



Element 


Distance from target 

[along nominal 
beam axis] 


^eff 


Max field/gradient 


Aperture 
Radius 

(Pole tip) 


1/2-aperture 

[H/V] 




(mm) 


(mm) 


(T) 


(T/m) 


(mm) 


mm 


mm 


Ql 
Q2 

Q3 
Dl 


3000.0 
4400.0 
5800.0 
7979.1 


888 
888 
888 
1038 


1.6 


1.6 
1.6 
1.6 


101.5 
101.5 
101.5 






254.0 


76.0 


Decay Solenoid 


12210.7 


5000 


5.7 


- 


57.5 






Proton absorber 


14880 


Plastic sheets 15, 29, 49, 54 mm 


GVAl 


15050.0 


Scintillation counter (0.02 Xq) 


D2 


15808.1 


1038 


0.85 


- 




254.0 


76.0 


BPMl 


16992.0 


Scintillating fibres (0.005 Xq) 


Q4 
Q5 
Q6 


17661.6 
18821.6 
19981.6 


660 
660 
660 


- 


2.3 
2.3 
2.3 


176.0 
176.0 
176.0 






TOFO 


21088.0 


Segmented scintillator (0.12 Xq) 


Ckova 
Ckovb 


21251.5 
21910.9 


Aerogel threshold Cherenkov (0.019 Xq) 
Aerogel threshold Cherenkov (0.031 Xq) 


BPM2 


24293.7 


Scintillating fibres (0.005 Xq) 


Q7 
Q8 
Q9 


25293.7 
26453.7 
27613.7 


660 
660 
660 


- 


2.3 
2.3 
2.3 


176.0 
176.0 
176.0 






TOFl 


28793.1 


Segmented scintillator (0.12 Xq) 


T0F2 
KL 

Tag counters 


31198.1 
31323.1 
31423.1 


Segmented scintillator (0.12 Xq) 
Lead + scintillator (2.5 Xq) 
Scintillation bars (0.06 Xq) 



The MICE beam line works with both positively and negatively charged muons. The change- 
over is achieved by physically disconnecting and swapping over the magnet power-supply cables. 
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Figure 8. Upstream beam line. The target may be seen in the foreground of the photograph. The magnets 
that make up the upstream beam Hne (Ql-3 and Dl) He to the left of straight 7 of the ISIS proton synchrotron 
which can be seen in the centre of the figure. The luminosity monitor, surrounded by its white, borated- 
polythene shielding may be seen close to the right-hand edge of the photograph. 

A simple system for monitoring and recording the operating polarity of the two bending magnets 
(and by inference the other beam-line elements) has been deployed fESlI consisting of a pair of 
Honeywell LOHET-II Hall-effect sensors that are operating past saturation when placed in the 
dipole bores, even at minimal field strengths, and thus return one of the two well-defined voltages 
corresponding to the two possible polarities of each magnet. 

3.2.4 The proton absorber 

With the beam line set to transport positive particles, a large flux of protons is observed to exit 
the decay solenoid. The range of proton momenta selected by Dl is such that a large fraction of 
the protons reaching D2 will be accepted by Q4-6. As a result, unless action is taken to remove 
protons from the beam as it enters D2, an unacceptably large flux of protons would be transported 
to TOFO. In the momentum range of interest, the energy loss per unit length of material traversed 
is larger for protons than for pions. Borated-poly ethylene sheets placed in the beam as it emerges 
from the decay solenoid are used to remove the protons. The boron loading of the polyethylene 
serves to increase the absorption rate of neutrons. The sheets have thicknesses of 15 nrni, 29 mm, 
49 mm and 54 mm and are raised and lowered into the beam using a cable-pull mechanism. The 
absorbers can be operated together or independently. When all four absorbers are placed in the 
beam, the 147 mm of material is enough to stop protons with momentum up to ^ 500 MeV/c. 
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Figure 9. Downstream beam line. The ISIS synchrotron lies behind the curved wall in the background. 
The downstream quadrupoles (blue, hexagonal magnets) may be seen just downstream of the white concrete 
shielding. The time-of- flight counter TOFl is seen immediately downstream of the final quadrupole. TOFl 
is followed by the KL and the prototype EMR. 

The choice of absorber thickness depends on the required muon-beam momentum. In a beam 
line set to deliver a nominal muon momentum of 140MeV/c, Dl is set to select 360MeV/c parti- 
cles. In this case, a proton-absorber thickness of 44 mm is required to reduce the proton contami- 
nation in the muon beam delivered to MICE to less than 1%. The beam-line tune that delivers the 
nominal muon momentum of 240MeV/c requires that Dl be set to select particles with a momen- 
tum of 507MeV/c, and the full 147 mm proton-absorber thickness is then required to reduce the 
proton contamination to the few-% level. Pure positive muon beams of intermediate momentum 
are delivered using appropriate combinations of the four sheets. 

3.2.5 Diffuser 

The final element of the downstream section of the beam line is the diffuser. By introduction of 
material with a large radiation length into the beam, multiple Coulomb scattering increases the di- 
vergence of the beam, enlarging the emittance and providing an approximately matched beam in the 
upstream spectrometer solenoid. Since the diffuser must operate in the field of the 4 T spectrome- 
ter, electric motors and other magnetic components cannot be used. The diffuser design consists of 
a stainless steel drum, which is inserted into the upstream section of the first spectrometer solenoid 
fE9l1 . The drum contains four irises, two of brass (2.97 mm, 0.2 Xq, and 5.94 mm, 0.4 Xq thick) 
and two of tungsten (2.80 nrni, 0.8 Xq, and 5.60 mm, 1.6 Xq). By selecting different combinations 
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of irises, one can achieve a total of 3Xo in steps of 0.2 Xq. Each iris consists of two planes of 
four segments of material, supported by a Tufnol fBDII ring of radius 100 mm (see figure [T^) and is 
operated by a non-magnetic, air-driven actuator. The position of each iris is monitored by a set of 
optical sensors. When closed, the iris inserts solid petals into the beam; when it is open, the petals 
are stowed away within the Tufnol ring. 



Irises 



Optical 
sensors 





Figure 10. Left panel: schematic of the MICE diffuser showing the irises. Right panel: one of the brass 
irises being tested. 



3.2.6 Beam-line diagnostics 

The downstream beam line is instrumented with a scintillation counter (GVAl), placed just after 
the proton absorber, upstream of D2. GVAl consists of a single slab of scintillator with 18 x 
18 cm^ active area and 1 cm thickness, read out by a 2 inch EMI photomultiplier. Two Beam Profile 
Monitors mounted on the upstream faces of Q4 (BPMl) and Q7 (BPM2) were used to measure 
particle rate and beam position. The BPMs each comprise two planes of Kuraray 1 mm scintillating 
fibres. BPMl has an area of 20 x 20 cm^, while BPM2 has an area of 45 x 45 cm^. The layers of 
optical fibres are supported by a sheet of Rohacell [Mil foam 2.55 mm thick. An aluminium box 
encases the BPMs, except in the acceptance region of the fibres where four layers of Tedlar fB^ 
keep the BPM light-tight. The light from the fibres is read out by two 64-channel multi-anode 
Burle fE311 8501 1-501 photomultiplier tubes (one per plane). Further details on both GVAl and the 
BPMs are found in table 0. 

3.3 Optimisation of the beam line 

The upstream and downstream sections of the beam line are loosely coupled, allowing each section 
to be optimised independently. In the upstream section, where pions travel through an evacuated 
beam pipe, the matrix-evolution code TRANSPORT [El, ESD was therefore used to perform the 
optical design. In the decay solenoid, pion decay must be simulated and in the downstream beam 
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Figure 11. MICE beam-line envelope for the baseline case with % = 67r mm • rad and p^ = 200 MeV/c. 
Left panel: pion optics for the upstream section. Right panel: muon optics for the downstream section. 
The lead scatterer ("diffuser") used in this simulation for practical reasons has been realized using brass and 
tungsten. For both panels the top half of the plot shows the vertical envelope and the bottom half shows the 
horizontal envelope. 



line, where there is no beam pipe, particle interactions in air and in the detectors must be taken into 
account. Therefore, the optics of the decay solenoid and the downstream beam line were developed 
using Decay-TURTLE [E2L E3]- The result of these simulations are illustrated in figure [TTI [EH]. The 
left-hand panel shows the 1 a beam envelope obtained using TRANSPORT for a 454 MeV/c pion 
beam generated at the target. The upstream quadrupoles (Ql-3) are configured in a focusing- 
defocusing-focusing (EDE) configuration in the horizontal plane (DED in the vertical plane). The 
particle rate is limited by the aperture of Q2 in the vertical plane and by the apertures of Ql and 
Q3 in the horizontal plane. The right-hand panel of figure |n| shows the 1 a envelope obtained 
using TURTLE to track a sample of muons (/7^ = 238MeV/c at D2) towards the diffuser. This 
muon beam was optimised to deliver muons that would have a nominal momentum of 200 MeV/c 
and a nominal emittance of 6 tt mm • rad. The optical designs were then further refined using two 
GEANT4-based ES, E3] simulation codes: G4beamline EO, SI] and G4MICE [H ED. 
Pion transmission in the upstream beam line was maximized by increasing the Ql-3 excitation of 
^ 10% with respect to that obtained from a TRANSPORT simulation. 



3.3.1 Downstream optimisation 

The beam optics are tuned for the experiment set up in the Step VI configuration. The downstream 
beam line is designed to provide an optically matched beam of desired nominal emittance ^^ and 
nominal momentum p^. The parameters available for tuning are: 

1. The thickness of the diffuser, which is used to tune the beam divergence in the upstream 
spectrometer solenoid; 

2. The quadrupole currents, which are used to tune the beam size and divergence, i.e., the 
incoming Twiss parameters (fif^^^ and a}^^ a}^) at the upstream face of the diffuser; 

3. The nominal momentum, which can be adjusted with the excitation of D2. 
Ideally, the conditions to be fulfilled simultaneously are: 
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1. The beam is matched in the upstream spectrometer solenoid, downstream of the diffuser. The 
matching conditions are: jS^ = jS^; = j8 = 2p2[GeV/c]/(0.35[T]) and a^^ 0Cy = a^O; 



The beam size at the entrance of the diffuser (essentially the same as the beam size at the exit 
of the diffuser) is 
with y^E/m^\ 



of the diffuser) is such that the rms width of the beam in x and 3; satisfies a^ ^Oy ^ {EnP ) //, 



3. The beam angular divergence at the upstream surface of the diffuser (oj, and Gy), increased 
by multiple Coulomb scattering in the diffuser (cr^j^^^), is equal to the desired angular di- 
vergence in the spectrometer solenoid, i.e.: O"^, + Gq ^^^^ = Oy + Gq y^^^ = £n/{P7)'^ 

4. The central momentum selected by D2 matches the nominal momentum p^, once corrected 
for the calculated energy loss in all material encountered by the beam as it travels from D2 
to the central absorber. 

The horizontal beam dispersion is neglected in the calculations. It is not possible to fulfil all 
conditions simultaneously because the emittance of the beam coming out of D2, and transported in 
the quadrupole channel, is different in the horizontal and vertical planes. Nevertheless, an average 
match can be obtained by using the 4D covariance matrix, the resulting 4D emittance and the Twiss 
parameters of the beam at the upstream face of the diffuser, and transporting the beam through the 
diffuser. That the beam is only approximately matched results in a somewhat elliptical beam in the 
spectrometer solenoid; however this affects the cooling performance only marginally. 

Beam settings for various combinations of (Sn^Pz) ^^e calculated for a matrix of values 8^ = 
(3,6, 10) TT mm • rad and p^ = (140,200,240) MeV/c. The calculation procedure \M\ starts by fix- 
ing the emittance, the Twiss parameters and p^ inside the spectrometer solenoid at the downstream 
face of the diffuser. The Twiss parameters and the momentum are back-extrapolated through the 
diffuser. The diffuser thickness, d, and the currents in quadrupoles Q4-6 and Q7-9 are varied to 
match the input beam, using G4beamline to simulate the MICE downstream beam line. The pro- 
cedure is first applied to the reference point (6, 200). This beam line is scaled with momentum to 
provide a first approximation, MO, of each of the other points in the emittance-momentum matrix. 
Then, each point is re-optimised to produce the final set of optimised beam lines. Ml. The values 
of the Twiss parameters, diffuser thickness and momentum at the upstream face of the diffuser for 
the nine points of the matrix are given in table 0. 

4. The particle identification system in MICE Step I 

With the exception of the beam-line diagnostics and the luminosity monitor described in sections 
|3.2.6| and |3.2.2| , the detectors installed along the beam line form part of the particle identifica- 
tion (PID) system of the MICE experiment. The PID system upstream of the first spectrometer 
solenoid is composed of two time-of-flight (TOP) stations (TOFO and TOPI) fHHl and two thresh- 
old Cherenkov counters (Ckova and Ckovb) fHTII . Together the two Cherenkov counters will pro- 
vide 7r//i separation up to 365 MeV/c. TOPO, Ckova and Ckovb are inside the Decay Solenoid 
Area (DSA). ^ 



^The DSA is a closed area within the MICE Hall, just outside the ISIS vault area, that contains the DS and the first 
particle identification detectors. 



-17- 



Table 4. (En^Pz) matrix for the MICE programme [H3I1 . The Twiss parameters (a, j8) are those required at 
the upstream face of the diffuser of thickness d for a matched beam. 







ft (MeV/c) 




140 


200 


240 






J = 0.0 


J = 0.0 


(/ = 0.0 




^ 


Prf;//=151MeV/c 


Pdiff=201 MeV/c 


p<;,y/=245 MeV/c 






a=0.2 


a=0.1 


a=0.1 






13=56 cm 


jS=36 cm 


J3=42 cm 




(i = 0.9Xo 


(i=1.3Xo 


d=l3Xo 


B 
B 
1^ 


6 


Paiff =156 MeV/c 


Pdiff=2l5 MeV/c 


Pdiff=256 MeV/c 


a=0.3 


a=0.2 


a=0.2 


i 




j3=113cm 


jS=78 cm 


j3=80 cm 




d=l.^Xo 


(i = 2.8Xo 


d = 2.8Xo 




10 


Pdiff =164 MeV/c 


Pdiff=229 MeV/c 


Pdiff=261 MeV/c 




a=0.6 


a=0.4 


a=0.3 






J3=198cm 


J3=131cm 


jS=129cm 



The TOP system is required to reject pions in the incoming muon beam with an efficiency 
in excess of 99%. In addition, the precision of the TOP time measurement must be sufficient to 
allow the phase at which the muon enters the RP cavities to be determined to 5°. To satisfy these 
requirements, the resolution of each TOP station must be ^ 50 ps. The two Cherenkov detectors 
have been designed to guarantee muon-identification purities better than 99.7% in the momentum 
range 210MeV/c to 365MeV/c fffSll . At lower momenta, 7r//i separation is obtained using the 
TOP measurement, both Cherenkov detectors being blind to both particle types. 

The identification of particles downstream of the second spectrometer solenoid is provided by 
a third TOP station (TOP2) fH9l1 and an electromagnetic calorimeter (EMC) by which muons may 
be distinguished from electrons (or positrons) produced in muon decays downstream of TOPI. The 
electromagnetic calorimeter consists of two parts: "KLOE Light" (KL), a lead- scintillator device 
based on the KLOE calorimeter design flT^ : and the "Electron-Muon Ranger" (EMR), a fully 
active, scintillator calorimeter QH]). The left-hand panel of figure |TJ shows the positions of TOF2 
and KL in the MICE muon beam line during Step I data-taking. The right-hand panel shows a 
CAD drawing of the downstream PID detectors (TOP2, KL, EMR) on their platform. The Ckov, 
TOPO, TOPI, TOP2 and KL detectors were installed in stages in the MICE Hall. The EMR was 
not available for Step I data-taking. 



4.1 The Cherenkov counters 

Two aerogel Cherenkov counters are employed fHTII . The muon threshold in the upstream counter. 



Ckova, is set at /^^^ = 278 MeV/c while the muon threshold for Ckovb is set at p^f^j^ = 210 MeV/c. 
Pions thresholds are p^^^ = 367 MeV/c and p^^^ = 277 MeV/c for Ckova and Ckovb respectively. 
The densities of the aerogels are pa = 0.225 gcm~^ and pt, = 0.370 gcm~^, with indices of re- 
fraction Ha = L07 and nt, = L12. In figure |T3| an exploded view of one counter is shown. The 
aerogel tiles are two layers thick (2.3 cm in total) and cover an area of 46 x 46 cm^. Pour 8 inch 
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Figure 12. T0F2 is shown in front of KL on the final downstream platform (left panel) and the downstream 
PID section including T0F2, KL and EMR on their platform (right panel). 

EMI 9356KB PMTs collect the Cherenkov light in each counter. Due to the high particle rate, the 
digitisation of the pulse profile is performed using a very high frequency sampling digitiser, CAEN 
V1731 [501 (1 GS/s maximum sampling rate). The sampling digitiser is connected directly to the 
PMTs through a coaxial cable. 

For the 140 MeV/c beams, both pions and muons are below threshold for both of the Cherenkov 
counters. For 200MeV/c beams, pions are below threshold for both Ckova and Ckovb, while 
muons are above threshold only for Ckovb. For 240MeV/c beams, pions are above threshold for 
Ckovb while muons are above threshold for both Ckova and Ckovb. Algorithms are being written 
to use the information from both counters to provide Tt^ji^e likelihoods. The likelihood information 
will be combined with the TOF measurement. Only the TOE system was used for PID in the Step 
I data analysis presented in the present paper. 

4.2 The TOF detectors 

All TOF stations share a common design based on 1 inch thick slabs of fast scintillator material. 
Bicron BC-420 scintillator El was used for TOFO and BC-404 was used for TOFl and TOF2. 
The slabs in each station form two planes with x and 3; orientations to increase measurement re- 
dundancy. Each slab is read out at each end by a Hamamatsu R4998 fast photomultiplier (rise time 
< 1 ns). The active areas of TOFO ,TOFl and TOF2 are 40 x 40 cm^, 42 x 42 cm^ and 60 x 60 cm^, 
respectively. The strip width is 4 cm for TOFO and 6 cm for the other two stations. All downstream 
PID detectors and the TOFl station will be shielded against stray magnetic fields of up to 1300 G 
(with a < 400 G longitudinal component) due to the presence of the spectrometer solenoids. TOFl 
will be shielded by a double- sided iron cage which fully contains the detector. The iron cage has 
a hole for the beam. The TOF2 and KL PMTs are shielded by soft iron boxes [22], as shown in 
Figure 0. In MICE Step I, only the local PMT shields for TOF2 and KL were installed. 

The TOF stations must sustain a high instantaneous particle rate (up to 1.5 MHz for TOFO). 
The rate capabilities of the R4998 PMT were tested in the laboratory with a dedicated setup based 
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Figure 13. Aerogel Cherenkov counter blowup: a) entrance window, b) mirror, c) aerogel mosaic, d) acetate 
window, e) GORE reflector panel, f) exit window and g) 8 inch PMT in iron shield. 




Figure 14. Shielding used for TOP PMTs. Left panel: TOFl global cage showing. A) extraction brackets, 
B/E) shielding plates, C) TOFl detector in working position, D) rails to move TOFl outside the shielding 
cage. Right panel: T0F2 local shielding with massive soft iron boxes. 

on a fast laser fHHl . The PMT rate capability was increased through the use of an active divider 
base. These signals, after a splitter, are sent to a leading-edge LeCroy 4415 discriminator followed 
by a CAEN VI 290 TDC for time measurement. The signal is also sent, after an RC shaper, to a 
sampling digitiser CAEN VI 724 EADC (lOOMS/s maximum sampling rate) in order to measure 
the pulse height for time-walk correction. 

The passive splitter is designed to match the impedance of the 50 Q. coaxial cable, coming from 
the PMTs, with the 120 f^ discriminator and shaper inputs. The shaping circuit is used to extend 
the duration of the short PMT pulse, so that it can be finely sampled by the digitiser. Software 
processing of the digitised pulse is needed for the charge measurement and time-walk correction. 

The signal arrival time, ti (/ = 1,2), at the photo-cathode of the left/right PMT of a scintillator 
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slab of a TOF detector plane is given by: 

,, = ,o + ^^^^^^ + 5,; (4.1) 

where ^o is the particle arrival time, x its distance from the counter centre, L the scintillator length, 
Veff the effective velocity of light in the scintillator slab and 5/ includes all time delays (in cables, 
PMT transit time, etc.). The transverse impact position, w, of a particle on a TOF station may be 
reconstructed from the difference between the time measurements from the two PMTs /, j at the 
ends of a counter as: 

u='^x{{u-5i)-{tj-5j)). (4.2) 

Transformation to the MICE coordinate system is straightforward: for vertical planes x = w and for 
horizontal planes y = u. The measured weighted average for Ve/f is 13.52 ± 0.30 cm/ns giving a 
spatial resolution cj^ = cjy = \/2 x Ve/f x Ot ^l cm, with Of the time resolution for the TOF station 
in question. 

The calculation of the delays 5/ (the "time calibration") is performed using Step I data as 
described in detail in fHHl and fl5311 . 

The use of leading-edge discriminators in the time-measurement electronics causes the threshold- 
crossing time to depend on the quantity of charge collected, an effect referred to as "time walk". 
The difference between the time measured by each TDC and a reference time was determined as a 
function of the maximum of the FADC signal and used to correct for time-walk. Pre-equalisation 
of the amplitude response of the scintillation counters in each TOF plane was performed so that the 
time-walk corrections for the two PMTs on a particular scintillator counter were similar. The pre- 
equalisation was performed using a YAP pulser unit from SCIONIX [5^ yielding ^ 20 counts/s 
from a YAP:CE scintillation crystal spot-activated with ^"^^ Am. Pulse-height spectra were recorded 
in the laboratory and used to derive appropriate high- voltage settings for the PMTs. 

The performance of the TOF stations was determined by measuring the difference between 
the time recorded by the horizontal (x) and vertical (y) planes, At^y, of each station. If it is assumed 
that the resolution of each plane is the same, the width of the distribution of At^y can be used to 
determine the time resolution of a single plane, <7^^^. Figure |T3| shows the At^y distribution for each 
TOF station. They all have remarkably good Gaussian shapes. The widths of the distributions are 
^ 100 ps, corresponding to TOF detector resolutions Ot ^ 50 ps. 

Figure |T^ shows distributions of the time-of- flight between TOFO and TOFl. The left panel 
represents data taken with a tt ^ /i beam. It has a small contamination of electrons and pions. 
Similar beams will be used to demonstrate ionization cooling. The right panel shows data taken 
with a calibration beam. In this beam configuration electrons, muons and pions fall into three well- 
defined peaks. Similar plots for the time-of-flight between TOFO and TOF2 are shown in figure |T7| . 

The reference muon beam is a standard setup of the beam-line magnets that has been used 
to monitor the stability of the data-acquisition system and the beam line. Each data-taking shift 
begins with a run that reproduces this reference setup. The reference muon beam runs can also 
be used to validate the consistency of the calibration procedure over the running period. Figure 
[TS| shows the variation of the TOF detector resolutions, as computed from the At^y distribution. 



-21- 




800 1000 

A t [ps] 

TOF2 resolution | 



TOF1 resolution! 


intrles 


14348 




Mean 


-6 76 


= 500 


- J 


RMS 


123.3 




M 


X^ / ndf 


128.8/95 


400 


'- / 1 


Constant 
Mean 


474.2 + 5.0 
-5.622 + 1.004 




'- j \ 


Sigma 


119.6 + 0.8 


300 








200 


~ \ 








100 


"l 1 1 1 1 1 1 1 1 hnHlrti/r 1 1 1 1 1 1 1 1 1 1 iX 


»lnnl.l 1 1 1 1 







ma -BOO -600 -4qo -200 



200 400 600 800 1000 

At [ps] 



Entries 

Mean 

RMS 



114S2 
47.6 
108.3 
98.42 ; 38 



Constant 429.1+5.1 
Mean 48.1+1.0 

Sigma 105.8+0.8 




200 400 600 800 1000 

A t [ps] 



Figure 15. Time difference Atjy between vertical and horizontal slabs in TOFO, TOFl and T0F2. The 
trigger is on TOFl. 
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Figure 16. Time of flight between TOFO and TOFl for a muon beam (left) and a "calibration" beam (right). 

for each reference run. The average time resolutions of TOFO, TOFl and TOF2 are 52.2 ± 0.9 ps, 
59.5 zb 0.7 ps and 52.7 ± 1.1 ps respectively fB311 . The slightly worse resolution of TOFl arises 
from some of the PMTs used on TOFl being of slightly poorer quality than the other tubes used 
in the TOF system. ^ The resolution of the TOFO station (4 cm wide slabs) and that of the TOF2 



^This feature was corrected later by refurbishing all PMTs of TOFl, obtaining a TOFl detector resolution comparable 
to those of TOFO and T0F2 Q. 
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Figure 17. Time of flight between TOFO and T0F2 for a muon beam (left) and a "calibration beam" (right). 
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Figure 18. Stability of the time resolution (in ps) of the TOF stations versus running time. The covered 
period is about one month of data- taking. 

station (6 cm wide slabs) are similar, showing that light path-length fluctuations are negligible. The 
stability of the TOF measurement for both upstream (TOFl-TOFO) and downstream PID (TOF2- 
TOFO) is of the order of ^ ±30 ps (see figure |T^ for details). This matches well the required 
resolution for PID (^ 100 ps). 

4.3 KL detector 

KL is a KLOE-type sampling calorimeter flTHl . composed of extruded Pb foils in which scintillating 
fibres are placed in a volume ratio Scintillator : Pb ^ 2 : 1 (see figure ^. Since the particle energy 
is lower in MICE than in KLOE, the ratio of lead to fibre has been reduced from that adopted in 
flTBl (^1:1), hence the name "KLOE Light". The fibres chosen are Bicron BCF-12 with 1 mm 
diameter, scintillating in the blue. The distance between two neighbouring fibres in the same layer 
is L35 mm. The distance between two layers is 0.98 mm, one layer being shifted by half the fibre 
pitch with respect to the next. The overall detector dimensions, including magnetic shielding and 
housing of photomultiplier tubes and voltage dividers, is approximately 120 x 4 x 160cm^. In 
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Figure 19. Stability of the time-of-flight of electrons between TOFO and TOFl (left) and TOFO and T0F2 
(right) versus run number. The covered period is about one month of data- taking. 
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Figure 20. Left panel: schematic layout of KL extruded lead layer and composite with fibres; Right panel: 
a photograph of a three-cell module of KL with Winston cone light-guides. 



figures and ^ a schematic view of one exploded KL module and the global layout of the KL 
assembly are shown. The active volume of 93 x 4 x 93 cm^ is divided vertically into seven modules, 
which are supported by an iron frame. The iron frame shields the PMTs from magnetic fields. KL 
has a thickness of 2.5 Xq and ^ 0.15 A/^^ From tests performed at the e^/e~ test beam facility 
BTF [5711 of INFN LNF, a time resolution of At ^ TOps/v^ and an electron-energy resolution, 
fully dominated by sampling fluctuations, of AE/E ^ 7%/y/E are obtained for electron energies 
between 75 MeV and 350 MeV. 

KL has 21 cells and 42 readout channels (one at each end of each cell). The light is collected 
by Hamamatsu R1355 PMTs with E2624-11 voltage dividers, providing differential output pulses 
on twisted pair cables with 120^1 impedance at 50 MHz. The signal from the PMTs is sent to a 
shaper module, which shapes and stretches the signal in time in order to match the sampling rate 
of the flash ADCs. The flash ADC modules are the same 14 bit CAEN V1724 used for the TOF 
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Figure 21. Exploded view of one KL module, showing the assembly of the fibre/lead module, light guides 
(A), mu-metal shielding (B), PMTs (C) and voltage dividers (D). 




Figure 22. Global layout of KL assembly. The exploded view shows the various custom made components 
for support and magnetic shielding: in yellow and dark blue the mechanical support and the PMT soft iron 
magnetic shields, in green the iron bars housing the PMT voltage dividers and in red the additional iron bars 
covering the Winston cones light-guides. 



stations. 
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In MICE StepI setup, KL is followed by three, 1 inch thick, 10 x 100 cm^ scintillator bars, 
placed vertically side by side behind the centre of the detector, to tag any particles that pass through 
KL. These "Tag counters" are used to measure the number of particles that pass through KL. They 
will be removed when the EMR calorimeter is installed. 

The TOP system was used to study the response of KL to different particle types. From the 
distance between the TOFl and TOF2 stations and the measured time-of- flight, the average mo- 
mentum between TOFl and TOF2 of the impinging particle was determined. For particles with 
momentum larger than 300MeV/c, TOF2-TOF1 time-of-flight has insufficient resolution to deter- 
mine the particle type, hence the TOFl-TOFO time-of-flight is used. To estimate the momentum 
at the front of KL, the energy lost in the material of TOF2 is taken into account. 

The procedure described above may be used for muons and pions. Since electrons and positrons 
are relativistic, it is not possible to use the time-of-flight method. The momentum of electrons and 
positrons is estimated using a model that takes into account all materials traversed by the electrons 
and the magnetic fields traversed. 

The data used to study the response of KL cover the entire range of momentum of interest in 
MICE. The ADC product is defined as the product of the digitised signals from the left and right 
sides of one slab divided by their sum: 

ADCprod = 2 X ADCieft X ADCright/(ADCieft + ADCnght) ; (4.3) 

where the factor of 2 is present for normalisation. The product of the two sides compensates for the 
effect of light attenuation. In the upper left panel of figure ^ the KL response to muons at various 
momenta is shown. The deposited energy reaches its minimum (^ 30MeV) at a momentum of 
300MeV/c (ionization minimum in lead). The upper-right panel shows the KL response to pions. 
In the plots, the abscissa represents the sum of the ADC product from all slabs in KL above a 
given threshold. The bottom panel of figure ^ shows the typical response of KL to electrons. The 
fraction of 80 MeV/c electrons which pass through KL is ^ 70%. 

The fractions of electrons, muons and pions passing through KL (the "KL transparency") are 
shown in figure 0. KL must act as a pre-sampler for the EMR, introducing a minimal perturba- 
tion to incoming muons and pions. The KL ADC product, plotted as a function of the time-of- 
flight for different particle types, is shown in figure ^. A clear separation between electrons and 
muons/pions is visible. 

Preliminary studies of KL response in different TOP time windows show a K contamination 
in the muon beam of the order of ^ 1%, after a suitable cut on the KL ADC product. These results 
will be the subject of a forthcoming paper. 

A raw time estimate may be made, using a simple linear interpolation of the rising edge of the 
sampling FADC signal. This provides some redundancy in the global time measurement. Prelimi- 
nary laboratory measurements have shown that a ^ 210 ps preliminary resolution may be obtained 
by selecting cosmic muons impinging at the centre of a KL cell. 

A study of the stability of KL during data-taking showed that there were no dead or noisy 
channels. The pedestals are stable with an rms of ^ 2.5 ADC counts, corresponding to ^ 0.2% of 
a MIP energy deposition. The KL response to muons and pions is stable in time (over a several- 
month period) to within ^1%. 
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Figure 23. KL response (normalised to the number of events) to muons for various incident momenta (top 
left panel), to pions (top right panel), with various momenta, and to 80 MeV/c electrons at the entrance of 
KL (bottom panel). For an explanation of "ADC product" see the text. 

5. The MICE Step I Trigger, Data Acquisition and Online/Offline Software 

5.1 Trigger and Data acquisition 

The Detector Data Acquisition (DDAQ) system of the MICE experiment must be able to acquire 
data at a rate of up to ^600 muons per spill. To fulfil this requirement, the Front-End Electronics 
(FEE) must digitise the analogue signal in less than 500 ns and store the digitised data in buffer 
memory. The time before the arrival of the next spill (^1 s) is used to read out the buffers and store 
the data. Data are transferred off-line to the remote mass storage unit for subsequent analysis using 
Grid computing technology fl5Hl . 

The acquisition of the data coming from the detectors is based on VME FEE interfaced to 
Linux PC processors via CAEN VME-PCI V2718 interfaces fBOl . These modules are connected 
via optical links to PCI cards embedded in the acquisition PCs and may sustain transfer rates up to 
70Mbytes/s. 

The software framework performing the management of all the independent readout processes 
and the event building, as well as the user interface, has been developed from the DATE package 
provided by the ALICE experiment at CERN [5911 . The event-builder machines receive sub-events 
from the various VME crates through the local network, put them together in memory and store 
them locally in a 3 TB, RAID6, disk buffer. These data are subsequently transferred to a second 
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Figure 24. Fractions of electrons, muons and pions passing through KL and reaching the Tag counters. 
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Figure 25. Time of flight, as measured from TOFO and TOFl, versus KL response as ADC product, using a 
300MeV/c "caHbration" beam, with trigger on TOFL 

local machine, from which they can be transferred asynchronously to the RAL Tier 1 centre for 
archival. As this transfer is carried out using Grid protocols, the data files are also immediately 
registered in a global file catalogue and replicated to other remote sites for off-line analysis. 

The Detector DAQ is strongly dependent on the Trigger System, which is divided into two 
parts. The readout can only be performed at the end of the spill, due to the time structure of the 
MICE beam. Digital information related to the passage of the particles through the MICE cooling 
channel during the spill must therefore be buffered in the FEE. The first part of the Trigger System 
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is responsible for the generation of the "Particle Trigger" signal. This signal triggers the digitisation 
of the analog signals received from the detectors. A coincidence of the PMTs connected to the same 
scintillation counter in one of the TOF stations generates the Particle Trigger. 

Logic units used to generate the Particle Trigger signals are LeCroy 4516 (CAM AC) logic 
modules fE311 . These use the output signal of the LeCroy 4415 discriminators, perform the needed 
logical AND operation between the two PMTs at the ends of the same counter and issue a logical 
OR output per group of 8 slabs. NIM logic is then used to produce a global OR signal for all the 
slabs of a TOF station. The resulting signal is distributed to all the front-end boards involved in the 
experiment. For the moment only TOFO or TOFl detectors may be used as a source of the Particle 
Trigger signal. 

The second part of the Trigger System generates the so called DAQ Trigger. This signal is 
generated after the extraction of the ISIS proton beam and causes the readout and storage of the 
digital data corresponding to all the particle triggers received during a spill. The DAQ Trigger is 
generated from the same signal as that sent to the target system to trigger the spill. 

5.2 The control system and online monitoring of MICE 

The EPICS (Experimental Physics and Industrial Control System) platform was chosen for the 
MICE controls and monitoring system because of its reliability, existing drivers for a wide variety 
of devices, flexibility and world-wide support network fEDII . The EPICS backbone is a local area 
network to which hardware components are interfaced via EPICS input/output controllers, while 
the user has access through EPICS channel access. 

All the equipment required for Step I of MICE, including the MICE Muon Beam, the TOF 
and Ckov detectors as well as the systems used for environmental monitoring of the MICE Hall, 
was included in the controls and monitoring system. The target, beam line, decay solenoid and the 
proton absorber and beam stop (which are manually operated) are monitored. 

The high- voltage system for the PID detectors is based on CAEN SY127 and CAEN SY527 
systems interfaced, via CAENET, with CAEN PCI A 1303 cards [1501 . The control system is used 
to set-up and monitor the high voltage and includes the functionality to turn the detectors on and 
off and set ramp rates, operating voltages and current limits. The environmental monitoring system 
includes the measurement of the temperatures for the TOF detectors and the temperatures and 
internal humidity of the two Ckov detectors. 

The MICE Hall is monitored using temperature and humidity probes, a neutron monitor and 
water-leak detectors. Selected parameters are used by the EPICS Alarm Handlers which compare 
read values to preset limits. The EPICS Archiver is used to archive selected parameter values, to 
be used later as input to correction algorithms being developed for the offline analysis. 

A MICE Configuration Database (ConfigDB) has been developed and will be used to store the 
parameters that determine the operation of the experiment as well as alarm limits for the various 
components EH, E3]. The TOF and Ckov systems read the parameter set- values and limits from 
the ConfigDB during initialisation. Alarms occur when parameter values drift beyond these limits. 

In addition to the Control and Monitoring system, the MICE reconstruction program runs 
online in parallel to data taking. This allows a detailed inspection of relevant physical parameters 
by the users on shift. 
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5.3 The MICE Offline Software System 

The simulation, design and analysis of the experiment are supported by a variety of software tools. 
Most useful are G4MICE [H ED and G4beamline [H, both of which make use of the GEANT4 
toolkit EBD- G4MICE provides an extensible framework for offline simulation, reconstruction and 
analysis as well as online monitoring and is the standard software in the experiment for all Step I 
analysis work. G4MICE contains configuration files that describe the geometry of the experiment 
(including detectors, magnets, absorbers and cavities), the cabling of the various detectors and their 
calibrations. Acting as a client to the Configuration Database, G4MICE is able to reconstruct data 
taken in any experimental configuration. 

In addition to the GEANT4 toolkit for high-energy physics simulation, G4MICE also exploits 
the Fermilab BeamTools package to describe magnets, absorbers and cavities. Inside G4MICE, the 
Class Library for High Energy Physics (CLHEP) package fl53l provides classes useful for random- 
number generation and other applications including mathematical functions (complemented by the 
GNU Scientific Library, GSL). In addition, the ROOT analysis framework fl5311 is used to develop 
graphical user interfaces and to perform data analysis using a wide variety of libraries including 
classes for drawing histograms, file persistency and graph fitting. 

G4beamline provides fast simulation of particle behaviour in a variety of electric and magnetic 
fields and materials for a given geometry and set of beam-line optics. G4beamline is implemented 
in C++, although its associated user interfaces are implemented in Java. Users interact with the 
application exclusively through sets of data cards describing the parameters and geometry of a run, 
and are not required to write their own C++ code. G4beamline also includes visualisation features 
which allow users to view and explore the geometries generated from their data cards, as well as 
to study the particle species and trajectories of a sample of particles at any point in space. The 
QGSP_BIC hadronic package of GEANT4, which adequately describes particle production in the 
MICE target, has been used by both G4MICE and G4beamline in the optimisation of the beam line. 

6. Characterisation of the MICE Muon Beam 

MICE took data in the Step I configuration in summer 2010. Data were collected in each of the 
nine (En^Pz) configurations defined in section 3 for both positive and negative beams. With the 
exception of a single broken PMT in TOFO, all detectors worked well. The MICE Muon Beam was 
operated at a maximum rate of 34 muons per target dip (in 1 ms) at a 1 V-ms beam loss. About 
10 million triggers were collected over the various beam configurations. The data were used for 
detector calibration, to assess the particle identification capabilities of each of the detectors and 
to characterise the beam line. The main results on muon rates and comparisons between data and 
Monte Carlo are summarised in the next section. A first measurement of the beam emittance, using 
the TOE detectors only, was also obtained fE^ and will be the topic of a separate publication fl5511 . 

6.1 Particle rates and beam profiles in the MICE beam line 

A dedicated study of particle rate in the MICE Muon Beam as a function of ISIS beam loss was 
performed. The particle rates observed in the GVAl, BPM2, TOFO and TOFl detectors were 
recorded along with the integrated SP7 beam loss. A series of runs, consisting of approximately 
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Figure 26. Total particle rates per spill vs induced ISIS beam loss for a negative n ^ ji beam, with a 3.2 ms 
spill gate (left), and for a positive n ^ ji beam, with a 1 ms spill gate (right). 



400 target dips each, were taken, each run corresponding to a particular target setting. The ISIS 
beam conditions were stable across a run. The mean beam loss and the mean particle rates were 
calculated by averaging the beam loss and particle rates recorded over the 400-dip sample. Figure 
P5| shows the dependence of the BPM2, TOFO and TOFl average particle rates as a function of 
the average beam loss. Negative k ^ ji beam-line optics with a 3.2 ms spill gate are shown in 
the left-hand panel of Figure ^ while positive n ^ }! beam-line optics with a 1 ms spill gate are 
shown in the right-hand panel. Straight-line fits are also shown and give an excellent description of 
the data S, [El. 

Muons were selected by the requirement that the time-of- flight between TOFO and TOFl be in 
the range 26.2 < A/ < 32 ns. Pion contamination of the muon sample was estimated using G4MICE 
giving ^ 1% for the negative beam and values between 3 to 5% for the positive one. In the future, 
the K contamination will be further reduced using the two Cherenkov counters (see section 4). 

Table ^ summarises the muon rate determined for each of the nine beam configurations. Errors 
are a quadratic combination of the statistical error of the fits and the systematic uncertainties arising 
from the event selection. The results presented in table ^ indicate a muon rate of ^ 5 /i/s per Vms 
for a 3.2 ms gate for the negative beam and between 17 and 34 /i/s per Vms for a 1 ms gate for the 
positive beam. Studies are underway to determine the maximum beam loss induced by the MICE 
target which is acceptable during normal ISIS operations. In parallel, various proposals are being 
considered by which the muon rate per Vms of beam loss can be increased, including adding a 
beam "bump" to move the ISIS beam towards the MICE target, increasing the acceleration of the 
target shaft, and re-optimising the target geometry or material. 

A comparison between the beam profiles observed at TOFO and TOFl and those predicted by 
G4beamline has been carried out in order to validate the MICE simulations. Small misalignments 
between the TOE detectors and the magnetic axis formed by the quadrupole triplets were observed. 
Figures ^ and ^ show a comparison of particle rates per TOE slab in the x plane (left) and the 3; 
plane (right) after applying the misalignment corrections. The figures show that the beam profile at 
TOFO is asymmetric in the x projection. Furthermore, the x-distributions are observed to be about 
20% wider in the data than in the simulation. These distributions reflect the pion momentum 
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Table 5. Muon track rates for the two polarities of the MICE beam Hne. Counts are normaUsed to the 
measured V.ms beam loss used to characterise the target depth. 
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Figure 27. Slab by slab comparison between simulation and data for the x (left-hand panel) and y (right-hand 
panel) planes of TOFO for the (6,200) matrix element, with misalignment corrections applied. 
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Figure 28. Slab by slab comparison between simulation and data for the x (left-hand panel) and y (right-hand 
panel) planes of TOFl for the (6,200) matrix element, with misalignment corrections applied. 



spectrum as produced in 800 MeV proton Ti collisions, which is unlikely to be exactly modeled in 
GEANT4. 
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7. Conclusions 

The MICE Muon Beam on ISIS at the STFC Rutherford Appleton Laboratory has been described. 
In addition, the design and performance of the MICE beam-line instrumentation has been pre- 
sented. The MICE Muon Beam has been shown to deliver muon beams over the range of momen- 
tum and emittance required to carry out the MICE programme. Approximately 30 muons per 1 ms 
spill (at a rate of ^1 Hz, for a 1 V- ms beam loss) for a positive muon beam has been achieved, for 
a negligible activation of the ISIS synchrotron. Preliminary results on the pion contamination of 
the beam shows values at the per-cent level. The MICE Muon Beam and the instrumentation will 
serve the experiment as it is built up in a number of Steps, making it possible to prove the principle 
of the ionization-cooling technique that is essential to the realisation of the Neutrino Factory and 
the Muon Collider. 
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